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ABSTRACT 

We report on the results of a radio interferometric observation of NGC 7027 
in the CO J = 2-1 and ^'^CO J = 2-1 lines. The results are analyzed with 
morpho-kinematic models developed from the software tool Shape. Our goal is 
to reveal the morpho-kinematic properties of the central region of the nebula, 
and to explore the nature of unseen high- velocity jets that may have created the 
characteristic structure of the central region consisting of molecular and ionized 
components. A simple ellipsoidal shell model explains the intensity distribution 
around the systemic velocity, but the high velocity features deviate from the 
ellipsoidal model. Through the Shape automatic reconstruction model, we found 
a possible trail of a jet only in one direction, but no other possible holes were 
created by the passage of a jet. 

Subject headings: stars: Planetary Nebulae — stars: carbon — stars: imaging 
— stars: individual (NGC 7027) — stars: kinematics — stars: winds, outflows 



1. Introduction 

Although many planetary nebulae (PNe) have simple apparent morphologies, their true 
3-D structures are in fact not well known. This is true even for the most well-observed 
PN NGC 7027. In the optical and radio images, it has a clear shape of a ring, similar 
in appearance to the Ring Nebula NGC 6720. Its optical and radio continuum structures. 
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corresponding to the ionized gas component of the nebula, have been analyzed by various 
authors, generally suggest i ng that it is an ellipsoidal shell inclined at ~ 30° to the line of sight 



flScottlll973l : lMassonlll989l : iRoelfsema et al.lll99ll : IWalsh fc Clegdll994l ). Observ ations in the 



2000 


). There is 


(Cox et al. 


1997) 



line, however, reveal that it has a bipolar structure ( iLatter et al. 



NGC 7027, being a very young PN (dynamical age ~ 600 years; iMassonI Il989l ). still 
possesses the remnant of the molecular envelope from the asymptoti c giant branch (AGB) 
progenitor. Its CO envelope was first detected bv lMufson et al.l (I19751) and has been imaged 



by single-dish telescopes and interferometric te c hniques ( Massoii et al.l Il985l : iBieging et al. 



199ll : Ijaminet et al.lll99ll : ICraham et al.lll993bt iFong et al. 



20061 ). In addition to high spa- 



tial resolution, the molecular-line interferometric observations have the additional advantage 
that the high spectral resolution observations provide the additional dimension of kinematic 
structure. If interpreted by a proper analytical tool, interferometric data can provide pow- 
erful insights into the true 3-D structure of the nebula. 

In this paper, we report on the result of the millimeter-wave interferometric observation 
of NGC 7027 in the CO J = 2-1 and ^^CO J = 2-1 lines using the Submillimeter Arra}fl 
(SMA) . These observational results are analyzed with a morphokinematic modeling tool 
Shape (ISteffen &: LopezI |2006| ) . Our purpose is to reveal the structure of the central region 
consisting of molecular and ionized components and to find possible trails of jets that might 
have passed through the central region. 

The outline of this paper is as follows. In Section 2 we describe the details of the obser- 
vation and data reduction. In Section 3 we present the observational results. In Section 4 we 
present the procedure and results of the morphokinematic modeling with Shape. In Section 5 
we discuss the obtained results, and finally the present research is summarized in Section 6. 



2. Observation and Data Reduction 

The radio interferometric observation of NGC 7027 was made on September 20, 2004 
with the SMA. The SMA data set was downloaded from the Radio Telescope Data Center 
(RTDC) operated by the Harvard Smithsonian Center for Astrophysics. The data were 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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obtained under good atmospheric conditions with a zenith optical depth at 230 GHz of 
0.045-0.055. The two CO hnes (i.e., CO J = 2-1 and ^^CO J = 2-1 at 230.538000 GHz 
and 220.398681 GHz, respectively) were included in the observed frequency ranges of 229.4- 
231.3 GHz (USB) and 219.4-221.3 GHz (LSB). The array used 6 x 6 m elements in the 
compact array configuration. The baseline length ranged from 10 to 70 m. The field of view 
of a single antenna was about 54". The observation was interleaved every 20 minutes with 
nearby gain calibrators, 2202+422 and 2015+371, to track the phase variations over time. We 
calibrated the data using IDL-MIR, which is a data reduction software package developed 
by the SMA project. The absolute flux calibration was determined from observations of 
Uranus, and was approximately accurate to within 15%. We followed a standard calibration 
method. The final map has an accumulated on-source observing time of about 6 hours. The 
single-sideband system temperature ranged from 90 to 120 K, depending on the atmospheric 
conditions and airmass. The SMA correlator had a bandwidth of 2 GHz with a resolution 
of 0.812 MHz. The velocity resolution was roughly 1.06 km s~^ at the observed frequency. 
The phase center of the map was taken at R.A.= 21^07'^0.1589^ decl.= +42°14'10.169" 
(J2000). Image p rocessesing of the data was performed with the MIRIAD software package 



( ISault et al.lll995l ). The robust weighting gave a 2.74" x 1.37" CLEAN beam with a position 



angle of 75.9°. 

We detected the continuum emission by integrating over a 1.0 GHz range in USB (using 
emission free channels), and the flux density of the emission is measured to be 1.3 Jy. This 
value is s ignificantly smalle r than previous measurements at (or close to) 230 GHz (3-4 Jy; 
see, e.g., iHafez et al.l 120081 ). This is presumably because the extended (dust) component 
of the continuum emission is resolved out by interferometer. Since the main concern in 
the paper is the molecular component of the central region, of which the flux is completely 
recovered, we do not go into specifics about the continuum emission. 



3. Results 

The total flux profiles of the ^^CO J = 2-1 and ^^CO J = 2-1 lines are presented in 
Figure 1. The peak intensity of the ^^CO J = 2-1 and ^'^CO J = 2-1 lines are 49.1 Jy and 
5.3 Jy, respectively. The line widths (at the intensity level) of the ^^CO J = 2-1 and ^^CO 
J = 2-1 lines are 46 km s~^ and 62 km s~^, respectively. The velocity integrated intensities 
of the ^^CO J = 2-1 and ^^CO J = 2-1 lines are 1029.6 Jy km s"^ and 165.0 J y km s-\ re 



spect ively. Compared to a single-dish observation with the SMT 10 m telescope (jZhang et al. 



20081), 21 % and 81 % of the single-dish fluxes in the ^^CO J = 2-1 and ^^CO J = 2-1 lines 



have been recovered respectively in the present interferometric observations. Taking into 
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account the uncertainty in the flux measurements, we could consider that almost all of the 
single-dish flux of the ^^CO J = 2-1 line has been recovered, although a significant fraction 
of the ^^CO flux is lost. The majority of the missing flux of the ^^CQ emissio n must orig- 
inate in the extended spherical envelope ( iMasson et al.lll985l : iFong et al.ll2006l ). Given the 
low optical depth of the ^^CO line, the ^'^CO flux must primarily come from the dense, small 
region surrounding the central ionized region, and we recovered the ^^CO flux almost totally. 
For this reason, we believe that the ^^CO flux emitted from the molecular component of the 
central region, which is the main concern of this research, must be almost fully recovered 
also. Note that in this paper the "central region" means the region consisting of the central 
ionized region and surrounding molecular component that is traced by the ^'^CO emission. 

Figure 2 shows the total-flux intensity maps of the ^^CO J = 2 -1 and ^^CO J = 2-1 
lines, superimposed on the L-band (A = 18-20 cm) continuum image (IBains et al.ll2003l ). It 
is obvious that the ionized region is embedded inside the molecular emission region. The CO 
feature appears have an axial symmetric structure with the symmetry axis lying along the 
northwest to southeast direction (position angle~135°). The intensity of the CO emission 
is weaker in the northwest and southeast edges, and the emission region of the L-band 
continuum spills into this CO-faint regions. This appeara nce is consistent w ith previous 
observations in the CO J = 1-0 line (see, e.g.. Figure 11 in lBieging et al.lll99l[ ). The weak 
features seen in the corners of the upper panel of Figure 2 are uncleanable features, which 
are due presumably to an extended envelope. The ^^CO J = 2-1 feature is slightly larger 
than the ^'^CO J = 2-1 feature in the northwest to southeast direction. The angular size 
of the emission region of the ^^CO J = 2-1 line is roughly 26" (major axis: northwest to 
southeast) and 21" (minor axis: northeast to southwest) at a 3a level. This c orresponds to 



the h near sizes of 3.8 x 10^^ cm and 3.1 x 10^^ cm at a distance of 980 pc (jZijlstra et al. 



20081 ). Similarly, the angular size of the emission region of the ^^CO J = 2-1 line is roughly 



22" (major axis) and 18" (minor axis), corresponding to 3.2 x 10^^ cm and 2.6 x 10^^ cm at 
a distance of 980 pc. 

Figures 3 and 4 show velocity channel maps of the ^^CO J = 2-1 and ^^CO J = 2-1 
lines, respective ly. The emission regions seen in Figures 3 and 4 are very simi l ar to the ^^CO 
J = 1 - maps (jMasson et al.l Il985l : iBieging et al.l Il99ll : iJaminet et al.l Il99ll : iGraham et al. 



1993bl : iFong et al.l l2006l ). even though the spherical component is missing in the present 
12C0 J = 2-1 and ^^CO J = 2-1 maps. As in previous CO J = 1-0 maps, the ^^CO 
J = 2-1 and ^^CO J = 2-1 features exhibit a systematic motion: the northwest part of 
the feature is bright in low velocity channels, whereas the southeast part of the feature is 
bright in the high- velocity channels. The position-velocity diagram giving in Fig ure 5 is also 
consistent with the ^^CO J = 1-0 observations (see, e.g. 
been interpreted by an expanding ellipsoidal shell. 



Graham et al.lll993bl ). which has 
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4. Morpho-Kinematic Modeling with Shape 

Since the flux of the central region is almost totally recovered both in the ^^CO J = 2-1 
and ^^CO J = 2-1 lines, the present interfero metric data should be useful to derive the 
morpho-kinematic properties of the molecular component of the c entral region. We have 



constructed models using the Shape software originally developed by lSteffen fc LopezI (120061 ) 
for the analysis of optical imaging spectroscopy data of planetary nebulae. In this paper, 
we have applied the software for radio interferometric data. Although it does not include 
hydrodynamical or radiation transfer calculations. Shape is useful to create synthetic images 
at different velocity channels for comparison with actual observations. Shape builds a three- 
dimensional morpho-kinematic model from the distribution of moving particles. Therefore, 
the Shape modeling is potentially valid assuming "optically thin". Although the ^^CO line 
is often optically thick, several studies of NGC 7027 show that ^^CO J = 1-0 line is not 
very opt ically thick, i.e., r < 1, and the intensity is approximately proportional to column 



density (jPeguchi et al.l llQQOt IJaminet et al.lll99ll ). Additionally, we used the data of the 
^■^CO J = 2-1, which is usually optically thin, for the analyses in conjunction with the data 
of the ^^CO J = 1-0 line. Therefore, the assumption about the optical thickness in the 
Shape model does not seem to have harmful effects on the results. 

Although Shape does not directly handle the concept of "excitation", what Shape is 
actually doing is equivalent to assuming constant excitation (in the scheme of the radiation 
transfer calculation). Therefore, practically we are assuming, in the Shape modeling, that 
the structure of CO gas (of the central region) is fully traced by the present CO observation. 

Shape provides two different ways to construct a model: "interactive modeling" and 
"automatic reconstruction". In the interactive modeling, we begin by creating an initial 
structure and velocity fleld for the nebula. Synthetic p-v diagrams and channel maps are 
then produced from the model which are subsequently compared with the actual observa- 
tions. The geometry and velocity structure are then interactively modifled until satisfactory 
convergence is achieved. In the automatic reconstruction. Shape automatically determines 
the spatial distribution and velocity of particles by directly scanning the observational maps 
without any assumption about the nebular structure (note that, however, we need to assume 
the velocity law). 

In this work, we modeled the central region using both of the aforementioned methods. 
In the interactive modeling mode, we assume an ellipsoidal shell as the structur e of the central 



region , as in previous work to interpret the ^^CO J = 1-0 maps (see, e.g., I Graham et al. 



1993bl ). We also use the automatic reconstruction to explore the possible three-dimensional 



morphology of the nebula. The details of each model are given in the subsections below. 
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4.1. Ellipsoidal Shell Model with Shape Interactive Modeling 



The values of the parameters of the ellipsoidal shell (inclination angles, major- and 
minor-axis radii, and systematic velocity) are summarized in Table 1. For the kinematic 
structure, we assume a constant velocity field , as is the c ase fo r the kinematic model of the 
spectra-imaging data of the Br7 emission by I Cox et al.l (|2002| ). Th e expansion ve l ocity w e 
adopt (22 km s~^) is higher tha n the value (1 4 km s~ ^) given by iGraham et al.l (jl993bl ). 
The expansion velocity given by I Graham et al.l (jl993bl ) was not fast enough in the present 
case, and therefore we increased it to 22 km s~^, which fits the present observation better. 
A p ossible reason f o r the d iscrepancy of the expansion velocities between the present work 
and I Graham et al.l (jl993bl ) is that the present CO J = 2-1 observation is more sensitive 
for high velocity regions than previous CO J = 1-0 observations. The expansion velocity, 
22 km s~^, is consistent with the line- width of the CO J = 2-1 line. The volume of the 
ellipsoidal shell is sampled with 20,000 particles, which creates a sufficiently smooth map for 
comparison with the observation. A circular beam with a diameter of 2" is used to convolve 
the modeled image to simulate the observations. The velocity channel maps of the model is 
given in Figure 6. 

This simple ellipsoidal model, of course, cannot perfectly reproduce the observations. 
However, through this simple modeling we can still see that the observed features near 
the systematic velocity can be well approximated by an expanding ellipsoidal shell. The 
deviation of the observed high-velocity components from the ellipsoidal shell model may be 
related to the possible existence of high- velocity flows (discussed later in Section 5). 



4.2. 3D Model with Shape Automatic Reconstruction 

In addition to the ellipsoidal-shell model, a full 3-D model using the Shape automatic 
reconstruction is also generated. The particle distribution of this model is derived from 
direct scanning of the observed position-velocity diagrams and channel velocity maps. For 
the velocity field, we assume a linear velocity law v{r) = 2. Or" where v is the expanding 
velocity and r is the distance from the central star. We determined the proportionality 
constant (i.e., 2.0) so as to produce the most symmetric 3-D structure. We found that a 
constant velocity (as assumed in the ellipsoidal-shell model) is unable to fit as well as a 
linear-velocity law. Although one might think that the Hubble-type velocity law suggests 
linear, continuous acceleration of the envelope motion with a distance from the central star, 
the velocity law can be also explained by ballistic mo vements, due to a sudden acceleration 



in the past plus free expansion since then (see, e.g., iBujarraball l2008l ). Volumes of 50,000 



particles each are needed to create smooth maps for comparisons with the ^^CO and ^^CO 
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observations. A circular beam with a diameter of 2" is used for convolution with the model 
images. 

The channel-velocity maps of the automatic reconstruction models for the ^^CO J = 
2-1 (green) and ^^CO J = 2-1 (orange) lines are presented in Figure 7. The velocity- 
channel maps produced from the automatic reconstruction are almost exactly the same as 
the observational map; this is not surprising because the model is constructed by scanning 
directly the observed channel maps. The close match between the channel maps of the model 
and the observations tells us that the model is properly created. We note that the ^^CO 
data used to create the automatic reconstruction is averaged over a 4.4 km s~^ range, which 
is slightly larger than the integration range used for the ^^CO data (2.0 km s~^). Therefore, 
the orange feature (^^CO) in Figure 7 is slightly more enhanced in comparison with the green 
feature (^^CO). 

The ^^CO J = 2-1 and ^^CO J = 2-1 lines exhibit different intensity distributions in 
Figure 7. For example, in the velocity range from 16 km s~^ to 28 km s~^, the emission region 
of the ^^CO J = 2-1 line is surrounded by the ^^CO J = 2-1 emission. These differences are 
presumably due to the difference of the optical depth of the two lines. If this is the case, we 
can consider that the two lines trace the structure of the CO gas component of the central 
region in a mutually complementary manner. 

Figure 8 shows the three-dimensional visualization of the CO-gas component of the 
central region in the 24 different directions. As in Figure 7, both the ^^CO (green) and ^^CO 
(orange) data are combined to create this diagram. The (0, 0) image is the one as seen along 
the line of sight. According to Figure 8, the structure of the CO gas component of the central 
region seems to be a patchy torus, and no high-velocity jets are found. The inner and outer 
radii and thickness of the torus are ~1.8, 6.2 and 4.4", respectively. When translated into 
linear sizes, these angular sizes correspond to 2.6 x 10^^ cm, 9.1 x 10^^ cm and 6.5 x 10^^ cm 
at 980 pc. The apparent symmetry axis of the torus is almost the same as that found in 
the ellipsoidal-shell model (see. Table 1). In the edge-on views of the torus [for example, 
(9, $) = (45°, —90°), (45°, —45°)], we also find non-negligible amount of CO gas located 
in the polar regions. The CO structure may therefore be interpreted as an ellipsoidal shell 
with multiple holes in the polar regio n rather than as a patchy torus. This CO structure 



resembles the H2 structure derived in (jCox et al.ll2002l ). 



We note that the 2;-direction of the model can be freely squeezed or stretched without 
altering the intensity distribution [the z-direction is the direction perpendicular to the page 
in the (0, 0) image in Figure 8]. The reason why we have assumed v = 2r" here is that 
the proportionality factor of 2 gives the most symmetric distribution of particles in the 
three-dimensional space. 
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5. Structure of the Central Region and the Existence of High- Velocity Jets 



High- velocity outflows are common in PNe. They are responsible for the creation of the 
high-density shells seen in PNe, and the interaction between these shells and the fast wind 
lead to high-tempera ture shocked gas tha t lead to thermal X-ray emission. The d etection of 



20011) and BD +30°3639 flKastner et al. 



X-ray s in NGC 7027 flKastner et al.ll 19941) as we U as other PNe such as NGC 6543 fIChu et al. 



20001 ) has been attributed to this mechanism. In 



the latter two ex amples, their respective fast winds have been detected directly through 
UV spectroscopy (ILeuenhagen et al.lll996l : iPerinotto et al.lll989l ). No fast outflow has been 
directly detected in NGC 7027, probably because of its high dust extinction. Shocks have also 
been suggested to explain the o verabundance of th e N2H"'" molecule and the underabundance 
of CS and HNC in NGC 7027 Jzhang et al.lboosi ). 



The existence of multipolar nebulae have led to the suggestion that the fast outflows 
in PNe may not be spherically symm etric but are highly coUimated in certain directions 



The structure of the molecular shell (IGraham et al.l Il993al : iKastner et al.l Il994j : I Cox et al. 



19971 ) around the ionized ga s in NGC 7027 has been suggested as passages of multiple 
collimated high- velocity jets (jCox et al.l 120021 ). If we assume that the jets create "holes" 
in the central region, we may be able to retrace the trail of the unseen jets based on the 
structural inform ation of the central region. This method was flrst applied for NGC 7027 by 
Cox et al.l (I2OO2I ). who suggested position-angles of three possible jets, based on the structural 
information of the H2 emission region. Here, we revisit this discussion by adding the present 
CO data to the previous H2 and Br7 data. Since NGC 7027 is an ionization-bounded nebula 



(see, e.g.. iRodriguez et al.ll2009l ). the emission regions of the Br7, H2 and CO lines construct 
a stratifled structure. Therefore, by adding the new information obtained by the present CO 
observation, we can get better structural information of the central region. 

In order to compare the three datasets of the Br7, H2 and CO lines, we constructed 
Shape automatic reconstruction models for the CFHT data of the H2 1-0 S(l) and Br7 lines, 
in the same manner as the present CO data. The channel ma ps of the H2 1-0 S(l) and 
Br7 lines were taken from Figures 3 and 6 in I Cox et al.l ( 12002| ). We assumed a constant 
expanding velocity of 20 km s~^ for the Br7 model and a linear ve l ocity- law {V cx r) for the 
H2 1-0 S(l) model following the assumptions made by ICox et al.l ( 120021 ) in their kinematic 
modeling. Since the prop ortional constan t for the linear velocity law for the H2 1-0 S(l) 
model was not specifled in I Cox et al.l ( 120021 ). we found an appropriate constant [f (r) = 2.4r"] 
through tr i al and error, so that the model could reproduce the feature seen in Figur e 11 in 
Cox et al.l (120021 ). In fact, these velocity laws are not inconsistent with IWalsh et al.l ( 1l997l . 
at the inner region (roughly within 3-4"); since Walsh's law is constructed on the basis of 
the motion of the central ionized region, it is natural that the difference is large at the outer 
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neutral region). 

In Figure 9, we present the three-dimensional representations of the Shape model com- 
bining the Br7, H2 and CO data, and found that the emission regions of the three different 
lines occupy the spatial positions in a mutually complementary manner, creating a single, 
patchy ellipsoidal shape (the individual structure of the H2 and Bij emission is given in 
Figures 12 and 13). 

Using the combined Shape models, we tried to find "trails" left behind by the passage of 
jets. Specifically, we inspected the Shape particle distributions in various different directions 
using the Shape 3D module, and tried to find directions in which the particle distribution has 
a doughnut-like feature with a hole at the center. In this p rocedure, we fixe d the positions 



angles at PA=4°, —28° and —53°, which have been given by lCox et al.l (|2002[ ). and adjusted 



only the inclina tion angle ( i.e., i nclination of the axis to the line of sight). The position 



angles given by ICox et al.l (120021 ) are presented in Figure 10, superimposed on the (0, 0) 
image in Figure 9. We found that the combined model of the three emission lines exhibits 
a hole only in the direction of (PA, Incl.) = (— 53°, 145°) [hereafter, direction A (red bicone 
in Figures 11, 12 and 13)]. Except for this direction, emission regions of the three emission 
lines mutually close the angles, and therefore no other holes at other angles were found. 

We then removed the Br7 component from the model, and repeated a same inspection 
to find a doughnut-like feature suggesting the trail of jets. The reason we removed the Brj 
emission from the model is that the CO and H2 emissions exhibit more patchy structure as 
compared with the Br7 structure (see. Figures 11, 12 and 13), and also that the CO and 
H2 emission trace the neutral region while the Br7 emission traces the ionized region. With 
the model of just the CO and H2 emission, we found two other viewing angles showing a 
doughnut-like feature, namely (PA, Inch) = (4°, 125°) [hereafter, direction B (yellow bicone 
in Figures 11, 12 and 13)] and (PA, Incl.) = (— 28°, 65°) [hereafter, direction C (blue bicone 
in Figures 11, 12 and 13)]. Particle distributions seen in these directions are presented in 
Figures 11, 12 and 13. Here, we note that the extended particle distribution of the Br7 line 
(see. Figure 13) is due to the constant velocity law we assumed here. With the constant 
velocity law, the position of a particle is not uniquely determined in the z-direction, and 
therefore the particle distribution shows the weak, extended tail. 

In the top-right panels of Figures 11, 12 and 13, we can clearly see the doughnut-like 
distributions of the particles, while in the other two directions (see, middle-right and bottom- 
right panels in Figures 11, 12 and 13) the behavior of the Br7 emission is clearly different 
from that of the CO and H2 emission: the Br7 feature shows a hole only in the direction A, 
while the CO and H2 features show holes in all directions of A, B and C. 
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in direction A, all the three components of Br7, H2 and CO features show a hole and 
we can be relatively confident that there has been a passage of a physical jet along that 
direction. In the directions B and C, however, the Br7 emission region has no corresponding 
holes. 

If we assume a UV-radiation field with a bipolar shape, this might be alternative expla- 
nation of the hole-formation, because the UV radiation might create a hole on the structure 
of CO and H2 components by the effect of photodissociation. However, an important point is 
that the UV radiation does not seem to create a hole in the structure of Br7. In the present 
observation and modeling, we can clearly sec a hole in the Br7 structure (sec, Figsure 13) 
in the direction of PA=— 53°, and therefore the possibility of the UV radiation must be 
excluded, at least, in this direction. 

6. Summciry 

In this paper we report the results of an SMA observation of NGC 7027 in the CO 
J = 2-1 and ^^CO J = 2-1 lines. We also performed a morpho-kinematic analysis with 
Shape. The main results of this research are summarized below: 

1. The flux of the CO-gas component of the central region is almost fully recovered both 
in the CO J = 2-1 and ^^CO J = 2-1 lines, and the CO structure of the central region 
is apparent by combining the CO and ^^CO data. 

2. The simple ellipsoidal shell model explains the intensity distribution around the sys- 
temic velocity, but the high velocity features deviate from this ellipsoidal model. 

3. Through the Shape automatic reconstruction model, we found a possible trail of a jet 
along the direction of (PA, Incl.)=(-53°, 145°). 

4. Holes in other directions are seen in the CO and H2 structure, but these holes are 
closed by the Br7 emission. 

These observations suggest that the central region of NGC 7027 has a torus structure, 
with the polar directions being cleared away by a high-velocity outflow. As the ionization 
front propagates into the outer regions, NGC 7027 will develop into a bipolar nebula in 
emission lines as NGC 6302 and NGC 2346. There is a possibility that there exist high- 
velocity outflows in other directions, which may lead to the formation of a multipolar nebula. 
However, in such an interpretation, we should ensure consistency with the present result. 
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Results of this work suggest that molecular-hne interferometry observations when interpreted 
by a morpho-kinematic model can provide valuable insights into the true 3-D structure of 
PNe. 
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Fig. 2. — Total flux intensity maps in the ^^CO J = 2-1 {upper panel) and ^'^CO J = 2-1 
(lowe r panel] lines, sup erimposed on the L-band (A =18-20 cm) continuum image (gray 
scale; iBains et al.ll2003l ). In the upper panel, contour levels are —3.0, 3.0, 15.6, 28.1, 40.7, 
53.3, 65.9, 78.4 cr, and the 1 a level corresponds to 1.5 x 10~^ Jy beam~^. In the lower 
panel, contour levels are —3.0, 3.0, 5.9, 8.9, 11.8, 14.7, 17.6, 20.6, 23.5 a, and the 1 a level 
corresponds to 6.3 x 10^'^ Jy beam~^. The FWHM beam size is located in the bottom right 
corners. The origin of the coordinate corresponds to the phase center. 



- 15 - 




io ,e -to n^o^o"^ 10 -10 "lb 6^-10 to, p -io" 
RA offset (arcsec) 



Fig. 3. — Velocity channel maps of the SMA compact array data in the ^^CO J = 2-1 line. 
The velocity width of each channel is 2 km s~^ and the central velocity in km s^^ is located 
in the top left corner of each channel map. The contours levels are —3, 3, 14.71, 26.43, 
38.14, 49.86, 61.57, 73.29, 85 a, and the 1 a level corresponds to 8.04 x lO'^ Jy beam~^ 
The dashed contour corresponds to —3 a. The FWHM beam size is located in the bottom 
right corner of the last channel map. The origin of the coordinate corresponds to the phase 
center. 
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Fig. 4. — Velocity channel maps of the SMA compact array data in the ^^CO J = 2—1 line. 
The velocity width of each channel is 4.4 km s~^ and the central velocity in km s~^ is located 
in the top left corner of each channel map. The contours levels are —3.0, 3.0, 7.7, 12.4, 17.1, 
21.8, 26.5 and 31.2, and the 1 a level corresponds to 2.19 x 10^^ Jy beam~^. The dashed 
contour correspond to —3 a. The FWHM beam size is located in the bottom right corner of 
the last channel map. The origin of the coordinate corresponds to the phase center. 



-17- 




radial velocity {km/s) radial velocity {km/3) 



Fig. 5. — Position-velocity diagrams of the CO J = 2-1 (left panels) and ^'^CO J = 2-1 
(right panels) lines. The cuts passing the origin have position angles of 68° and 158°. The 
contours in the left panels are 0.34, 1.49, 2.64, 3.79, and 4.94 Jy beam~^. The contours in 
the right panels are 0.066, 0.116, 0.166, 0.216 and 0.266 Jy beani"^. 
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Fig. 6. — Channel maps of the best-fit eUipsoidal shell model. For simulating image convo- 
lution, a circular beam with a diameter of 2" is used. 
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Fig. 7. — Channel maps of the automatic reconstruction model. The green and orange 
colors respectively represent the ^^CO J = 2-1 and ^^CO J = 2-1 lines. For simulating 
image convolution, a circular beam with a diameter of 2" is used. 
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Fig. 8. — Three-dimensional representation of the molecular core of NGC 7027 in the ^^CO 
J = 2-1 and ^^CO J = 2-1 lines as seen from 24 directions separated by 45°. The line of 
view is identified by (O, $), where B is the zenith angle and $ is the azimuthal angle. The 
(0, 0) image is the nebula as seen for the Earth. The green and orange colors respectively 
represent the ^^CO J = 2-1 and ^'^CO J = 2-1 lines. The arrows indicate representative 
features. If the feature is lying in the far side and is obscured by another foreground feature, 
the position of the hidden feature is indicated by the dotted arrows. 
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Fig. 9. — Three-dimensional representation of the emission regions of the CO [^^CO J = 2-1 
and ^^CO J = 2-1], H2 [1-0 S(l)] and Br7 hnes. The notation of the hne of view is as in 
Figure 8. The CO, H2 and Br7 hnes are represented by the yeUow, blue and white colors, 
respectively (Note: both the ^^CO and ^^CO are represented by the same yellow color). 
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Fig. 10. — Position-angles of three possible jets as suggested by 
posed on the (0, 0) image in Figure 9. 
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Fig. 11. — Three-dimensional location of the three possible jets (yellow, blue and red bi- 
cones) seen from six different directions, embedded in the Shape model of the CO lines (pink 
particles). The CO data includes both the ^^CO and ^'^CO lines. Inclination angles (incli- 
nation to the axis to the line-of-sight) of the yellow, blue and red bicones are 125°, 65° and 
145°, respectively. Lengths and opening angles of the bicones are arbitrary. 




Fig. 12. — Similar diagrams to Figure 11 about the H2 1-0 S(l) line (blue particles). 




Fig. 13. — Similar diagrams to Figure 11 about the Br7 line (white particles). The elongated 
feature is not real (see, text). 
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Table 1. Ellipsoidal Shell Model of the Molecular Core 



Parameter Value 



Inclination of major axis to the line of sight 60° 

Position angle of major axis 155° 

Major axis radius 16" 

Minor axis radius 1" 

Thickness of the shell 3" 

System velocity (LSR) 25.6 km s"-*^ 

Expansion velocity 22 km s~^ 



